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Local pseudopotential (LPP) is an important component of the orbital free density functional 
theory (OF-DFT), which is a promising large scale simulation method that can still maintain in¬ 
formation of electron state in materials. Up to date, LPP is usually extracted from the solid state 
DFT calculations. It is unclear how to assess its transferability while applying to a much different 
chemical environment. Here we reveal a fundamental relation between the first principles norm- 
conserving PP (NCPP) and the LPP. Using the optimized effective potential method developed for 
exchange functional, we demonstrate that the LPP can be constructed optimally from the NCPP for 
a large number of elements. Our theory also reveals that the existence of an LPP is intrinsic to the 
elements, irrespective to the parameters used for the construction. Our method provides a unified 
method in constructing and assessing LPP in the framework of first principles pseudopotentials. 

PACS numbers: 71.15.Dx,71.15.Mb,31.10. +z,31.15.E- 


I. INTRODUCTION 


In recent years, orbital-free density functional the¬ 
ory (OF-DFT) 1 ’ 2 has attracted increasing interests due 
to its capability of simulating thousands to millions of 
atoms. 3-10 Comparing with other large scale simulation 
methods, OF-DFT maintains the electronic structure in¬ 
formation and is potentially applicable for systems un¬ 
dergo large chemical changes. However, OF-DFT has not 
become a mainstream method for large-scale simulation, 
due to the lack of both an accurate kinetic energy density 
functionals (KEDFs) and local pseudopotentials (LPPs) 
that are highly transferable in different chemical environ¬ 
ments. The recently developed nonlocal kinetic energy 
functional showed promising results for both metals and 
semiconductors. 11 ’ 12 On the other hand, it is still unclear 
how to proceed in systematically constructing transfer¬ 
able LPPs. Although several families of empirical LPPs 
are already available, 13-21 most of them are constructed 
to fit the solid state DFT results and only work for a 
small variation of chemical environment. 

One important concept in developing first prin¬ 
ciples pseudopotentials is the norm-conserving (NC) 
condition, 22,23 

[ nf CPP (r)d 3 r = [ nf E (r)d 3 r, (1) 

JQ J q 

in which n^ CPP (r) and n^ E (r) represent the orbital den¬ 
sities of NC pseudopotential and all-electron solutions, 
and i is the angular momentum number. This condi¬ 
tion requires the pseudo-charge enclosed within the core 
radius for each orbital must be identical to that of all¬ 
electron results and therefore ensures the transferabil¬ 
ity since it conserve the change of wavefunctions versus 


energy: 
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in which R(r, e) is the radial wavefunction of angular mo¬ 
mentum l. The thus constructed NCPPs are different 
for different angular momentum channels (orbitals) and 
therefore are non-local. Apparently, the NCPPs can not 
be used in OF-DFT calculations since the electrons in 
an orbital free scheme should feel the same potential. 
It is a long time question whether LPPs can be con¬ 
structed while still satisfying the NC condition. Here, 
we try to tackle this problem by altering it into the fol¬ 
lowing question: how can we optimally construct LPPs 
that are orbital independent from the NCPPs that are 
orbital dependent. 

We would like to point out the similarity between con¬ 
structing LPPs and obtaining the local exchange func¬ 
tional from its explicit Hartree-Fock form. As defined in 
Hartree-Forck method, the exact exchange energy is the 
functional of single particle wave functions, i.e. 

pexact 1 V [ ,3 * ^ (^ )</>* (^2 (r 2 )^ (u) 
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For systems without spin-orbit interaction, spin freedom 
is trivial to add and therefore is omitted throughout the 
paper. The corresponding exchange potential is non-local 
and orbital dependent. 


KMn) 



\n - r 2 1 


M r 2) 


( 4 ) 


Slater proposed a density average of Hartree-Fock nonlo¬ 
cal orbital dependent potentials which is known as Slater 





2 


potential. 24 
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turn channel. Correspondingly, the total ion-electron in¬ 
teraction energy can be written as 


( 5 ) 
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Tii(r)vi(r)d 3 r. 


( 8 ) 


Afterwards, a method that construct the optimal local 
exchange potential variationally from the explicit orbital 
expression of the exchange energy was established and 
applied for many cases. 25-27 


The expected local pseudopotential should reproduce the 
above total energy for different electron configurations of 
an atom except a constant shift. Such local potential can 
be calculated by taking a density derivative of Kj_ e , i.e. 
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This optimized effective potential (OEP) method can be 
transformed into an Algebraic equation, 


v(r) = 5Ei- e /6n{r). (9) 


Applying the OEP scheme and the KLI approximation 
originally developed for HF exchange potentials, one can 
obtain the corresponding KLI equation 25,26 as: 
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It can be simplihed to an algebraic Krieger-Li-Iafrate 
construction(KLI) 25 by neglecting the second term 
— V • ('0*(r)V^i(r)). The Slater potential is the 

first and the major term of it. It is worth to notice that 
OEP can be derived directly from the fact that the OEP 
orbitals are the first order shifts of the HF orbitals while 
conserving the charge density. 28,29 

We will show in this paper that the same idea of op¬ 
timally constructing local and orbital-independent ex¬ 
change potentials from the exact orbital dependent ex¬ 
change potential can be used to construct local pseudopo¬ 
tentials. The essential issue is the conservation of the 
NC condition, which as shown by our work, can be sat¬ 
isfied for large number of elements in the periodic table. 
We will also show that while NC condition is preserved, 
the Slater potential will be identical to KLI potential, 
which is an excellent approximation of the exact OEP of 
the NC pseudopotentials. On the other hand, our work 
also reveals that the transferable LPPs can not be con¬ 
structed for many elements, despite whichever scheme 
is used. Rather, as we demonstrated, the existence of 
a highly transferable LPP is an intrinsic property of an 
element. 


II. OPTIMIZED EFFECTIVE 
PSEUDOPOTENTIAL 


where the quantities with a bar above them are the 
orbital averages of the potentials, v(r) and Vi(r), i.e. 
v l = f rii(r)v(r)d 3 r, v\ = / rii{r)vi{r)d 3 r , and the sum¬ 
mation of i goes through all the angular momentum chan¬ 
nels. Again, the spin freedom is not considered and the 
spin index is omitted here. The first term is the long 
range Slater averaged potential, usiater and the second 
term is a short range correction. The derivation is similar 
to that of KLI equation for local exchange potential. 25,26 
It relies on the fact that the wavefunctions of local NC 
pseudopotential are the first order perturbation of the 
wavefunctions of semi-local NC pseudopotential. 

The above KLI construction will change the potential 
outside the core because of the non-Slater term. This 
is undesirable because the pseudopotentials outside the 
core should be identical to the true ionic potential. One 
possible solution is to construct KLI potential only inside 
the core and keep the original ionic potential outside the 
core. However, the resulted potential will be discontinu¬ 
ous at the core boundary. On the other hand, it is not 
hard to notice that the Slater potential that is the first 
and the major term of KLI potential does not change 
the potential outside the core region. This is because the 
first principles NC pseudopotentials outside the core are 
identical for different channels, and the Slater average is 
a density average of all angular momentum channel. 

The further consideration to the above problem relies 
on the relation between the NC condition and the OEP 
and KLI equations, which will be examined here. De¬ 
noting the orbital densities calculated from local pseu¬ 
dopotential and first principles NC pseudopotential as 
n 4jPP (r) and n^ CPP (r), the preservation of NC condition 
requires 

f n ppp {r)d 3 r = [ nf CPP (r)d 3 r, (11) 

Jn Jn 


We assume first principles NC pseudopotentials Vi(r ) 
for an atom, which are different to each angular momen- 


where Q denotes the space inside the core radius. Con¬ 
sidering that the change of the potential will make a first 
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TABLE I: Constructing and testing OEPP of Ga with various atomic configurations, eigenvalues of pseudo-atom and core 
radius, and with and without nonlinear core correction. The defininition of p a , p p and 8 P can be found in Appendix B 


Methods 

Configuration 

Rnlcc (cL-U) 

Rcut (s Ip) ( a.u ) 

e4eV) 

e P (eU) 

EZieV) 

Ps 

c 

Pp 

< 5 p 

TM-NCPP 

4s 2 4p 1 

- 

2.75 /2.75 

-9.1750 

-2.7384 

-58.0460 

0.7200 

0.3924 

0.0000 

OEPP 

4s 2 4p 1 

- 

2.75 /2.75 

-9.1966 

-2.7500 

-58.1033 

0.7250 

0.3807 

0.0182 

OEPP(<5 = 0.0018a.u) 

4s 2 4p 1 

- 

2.75 /2.75 

-9.1791 

-2.7414 

-58.0189 

0.7244 

0.3789 

0.0244 

TM-NCPP 

4s 2 4p 1 

1.75 

2.75 /2.75 

-9.1750 

-2.7384 

-104.8345 

0.7200 

0.3924 

0.0000 

OEPP 

4s 2 4p 4 

1.75 

2.75 /2.75 

-9.1944 

-2.7512 

-104.8916 

0.7248 

0.3808 

0.0172 

TM-NCPP 

4s u& 4p 2S 

- 

2.2 /2.2 

-10.7706 

-3.8364 

-47.9618 

0.5085 

0.2480 

0.0000 

OEPP 

4s°' 5 4p 2 ' 5 

- 

2.2 /2.2 

-10.8364 

-3.8418 

-48.0362 

0.5290 

0.2451 

0.0100 

TM-NCPP 

4s U0 4p 20 

1.75 

2.2 /2.2 

-10.7706 

-3.8364 

-93.3344 

0.5085 

0.2479 

0.0000 

OEPP 

4s°' 5 4p 2 ' 5 

1.75 

2.2 /2.2 

-10.8292 

-3.8421 

-93.4084 

0.5285 

0.2451 

0.0094 


TABLE II: s and p energy levels and the excitation energies for several atomic configurations of Mg, Ga and Sb, calculated by all 
electron potentials, Troullier-Martins pseudopotentials, OEPPS as well as the BLPS pseudopotentials. The two Troller-Martins 
pseudopotentials are constructed at the two configurations for OEPP and BLPS, respectively. Units are in eV. 





AE 

TMoepp 

TMblps 

OEPP 

BLPS 

Mg 



3s 2 3p u 

~~3s T 3p ! ~ 

3s 2 3p u 

3s 1 3p 1 

3s 2 3p u 


s 1 p 1 

s 

-5.7701 

-5.7702 

-5.7391 

-5.9200 

-5.8323 



V 

-2.1295 

-2.1295 

-2.1268 

-2.0873 

-2.1230 



A E 

3.5233 

3.5367 

3.5105 

3.6926 

3.5876 


s 2 p° 

s 

-4.7878 

-4.8104 

-4.7878 

-4.8216 

-4.8058 



V 

-1.3773 

-1.3760 

-1.3773 

-1.2668 

-1.3382 



A E 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


sY 

s 

-11.5278 

-11.4939 

-11.4416 

-11.5977 

-11.4980 



V 

-7.1642 

-7.1456 

-7.1272 

-6.9529 

-7.0234 



A E 

8.0742 

8.0824 

8.0459 

8.1341 

8.0802 

Ga 



4s 2 4p i 

4s 2 4p 1 

4s 2 4p i 

4s‘ ! 4p L 

4s 2 4p 1 




— 

r c = 2.75 

— 

8 = 0.0018 

— 


s 1 p 2 

s 

-10.2808 

-10.2781 

-10.2511 

-10.3780 

-10.0103 



V 

-3.5000 

-3.5027 

-3.5015 

-3.5368 

-3.4351 



A E 

6.6124 

6.6118 

6.5993 

6.6461 

6.4143 


s 2 p 1 

s 

-9.1750 

-9.1750 

-9.1750 

-9.1791 

-8.9113 



V 

-2.7384 

-2.7384 

-2.7384 

-2.7414 

-2.6713 



A E 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


s 1 p 1 

s 

-17.7538 

-17.7185 

-17.6771 

-17.8063 

-17.4160 



V 

-10.2007 

-10.1656 

-10.1639 

-10.1323 

-10.0634 



A E 

13.3385 

13.3279 

13.1192 

13.3720 

13.0427 

Sb 



5s 2 5p 3 

5s i 5p 4 

5s 2 5p a 

5s i 5p 4 

5s 2 5p a 


s 1 p 4 

s 

-13.8933 

-13.8933 

-13.8761 

-13.3325 

-12.2457 



V 

-5.5668 

-5.5668 

-5.5657 

-5.6422 

-5.4695 



A E 

8.2094 

8.2021 

8.2027 

7.5444 

6.6853 


s 2 p 3 

s 

-13.0893 

-13.0731 

-13.0893 

-12.4615 

-11.5606 



V 

-4.9991 

-4.9998 

-4.9991 

-5.0675 

-4.9692 



A E 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


sV 

s 

-21.7887 

-21.7689 

-21.7467 

-21.2557 

-19.9793 



V 

-12.8751 

-12.8494 

-12.8542 

-12.9483 

-12.6818 



A E 

17.3400 

17.3237 

17.3259 

16.7555 

15.6744 


order perturbation to the wavefunctions, one can prove 
the following relation for the potentials if the NC condi¬ 


tion Eqn.ll is satisfied (see Appendix A for proof), 

[ n PFP (r)(v(r) — Vi(r))d 3 r = 0. (12) 

Jn 

Eqn.12 can also be written as v l — v\ = 0. It suggests that 

















4 


if the NC condition is retained during the construction of 
local pseudopotential, the second term of the KLI poten¬ 
tial will become zero, indicating that the Slater potential 
is identical to KLI potential. 

KLI potential is an approximation of exact OEP poten¬ 
tial by neglecting a term involving the diversity of orbital 
change. 25,26 ^ (r)V<(>i(r)), in which = ipi~ 

i 

and ipi and (j)i are the wavefunctions of the original NC 
pseudopotentials and the constructed local pseudopoten¬ 
tials. This simplification can be interpreted as a mean- 
field approximation, since the neglected terms averaged 
over the ground state density p(r ) is zero, or equally the 
integration of the neglected term over space vanish, i. e. 

J=/d 3 r^V-(V ; *(r)V0 i ) = 0 (13) 

i 

If the construction of local pseudopotential retains the 
NC condition, this integral should be 0 for each individual 
orbital. 

h = f d\S7 • (#(r)V&(r)) 

Jq 

= ® dS r ij}*(s)V(/)i(s) ■ n = 0 (14) 

J n 

While the NC condition Eqn.ll is satisfied, the local 
pseudopotential will generate pseudo-wavefunctions that 
are identical to the wavefunctions of all electron poten¬ 
tials and the NC pseudopotentials outside the core re¬ 
gion, which means i/ii and should be exactly 0 at the 
core boundary and in the region out of the core. There¬ 
fore, we have shown that an optimized local pseudopo¬ 
tential can be constructed by taking the Slater average 
of the NC pseudopotentials if and only if the NC condi¬ 
tion is kept during the construction. If the NC condition 
can be retained, the local Slater averaged pseudopoten¬ 
tial is identical to the KLI potential of the semi-local NC 
pseudopotentials, which is very close to the exact OEP. 

III. CALCULATION DETAILS 

The FHI98 code 30 is modified to generate and test the 
proposed pseudopotentials. For comparison, the TM- 
NCPP 23 potentials are also generated by using FHI98 
code. The details of LPP construction for a set of 27 
elements, including the atomic configuration, the cutoff 
radius etc. are listed in Tabic SI. In order to improve 
the transferability, the nonlinear core-valence exchange- 
correlation 31 is included for some of the elements. The 
OEPP cutoff radii are adjusted to minimize the values of 
Si , S p and Se l (the definition of these quantities is given 
below ). The Kleinman-Bylander 32 form of pseudopoten¬ 
tials are used in solid calculations and we truncates the 
angular momentum at l max = 2. 

The structural relaxations and bulk property calcula¬ 
tions were carried out using KS-DFT as implemented 


in the CASTEP code. 33 Both local density approxi¬ 
mation (LDA) 34 exchange-correlation functionals and 
generalized gradient approximation (GGA) in Perdew- 
Burke-Ernzerhof (PBE) 35 form are used for exchange- 
correlation functionals. The appropriate energy cutoff 
and Monkhorst-Pack k meshes were chosen to ensure that 
enthalpy calculations for each system were well converged 
to less then 0.5meV/atom. A Fermi-Dirac smearing with 
a width of 0.1 eV is used for all metal systems. 

OF-DFT test calculations are carried out using our re¬ 
cently developed ATLAS software, 36 which is based on an 
efficient real space finite-difference method. We employ 
the Wang-Govind-Cater KEDF 3 ' with the parameters: 
7 = 2.7, « = (5 + V5)/6 and /3 = (5 - V5)/6. These 
parameters have been demonstrated to work well for Mg 
and A1 systems. 21,37 In finite-difference OF-DFT calcula¬ 
tions, the order of finite-difference ( N ) and the real space 
mesh gap h (determine the max plane-wave cutoff en¬ 
ergy E cut in the reciprocal space E cut = 7r 2 /2 h 2 ) should 
be tested for different systems. In our calculations for 
Mg and A1 systems, both energies converge to less than 
O.lmeV/atom, while setting N = 4 and h = 0.2A 


IV. TEST RESULTS 
A. Atom level 

In order to check whether a local NC pseudopoten¬ 
tial can be constructed for various elements, we define 
Si = f n m(r)(vi(r ) - vsiater(r))d 3 r as a measurement of 
the deviation from NC condition for the Slater averaged 
pseudopotential. According to Eqns. (11) and (12), 
if Si =0 then the NC condition is retained. Further¬ 
more, if Si is small, the changes of the orbitals are also 
small. Si estimates the energy changes for each orbital, 
i.e. £ 5PP — ef 101313 « Si, therefore they are good measures 
for the accuracy of LPPs. 38 

In the current work, we are going to focus on main 
group elements that only contain s and p electrons in 
their valence. For some p-block elements such as Ga and 
In, the d levels are close to the valence. Although they 
are completely filled, they may have effects to chemical 
bonds, and are included in the valence in many first prin¬ 
ciples pseudopotential constructions. However, the omis¬ 
sion of them in the valence usually will not yield errors 
that exceeds the use of OF-DFT and the local pseudopo¬ 
tentials. Therefore, we keep the d electrons in the core 
for all s and p elements. 

In the sp systems, we need to use only one parameter 
because <5o = — A (NC deviation value). Figure 1 shows 
<5o for elements from Li to Br. The Slater potentials are 
constructed from four kinds of NC pseudopotentials in¬ 
cluding BHS, 39 Kerker, 40 Vanderbilt 41 and TM. 23 In gen¬ 
eral, S 0 becomes smaller with increasing atomic number. 
However, for the atoms in the same row in periodic ta¬ 
ble, the higher the atomic number, the smaller the <5 0 
value. For most of the elements, <5 q does not depend on 
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FIG. 1: (Color figures available online.) The NC deviation val¬ 
ues <5o for elements from Li to Br. The original pseudopotentials 
are constructed by Troullier-Martins (TM), Kerker,Vanderbilt and 
BHS methods, and are represented by open squares, green trian¬ 
gles, red crosses and blue circles, respectively. 


the type of pseudopotentials except the noble gas atoms. 
This indicates the capability of constructing LPP is in¬ 
trinsic to elements and do not depend on the method of 
construction. 

Although OEPP is constructed directly from NCPP 
through the Slater average, we still have freedoms while 
constructing the NCPPs, which we can utilize to opti¬ 
mize the performance of OEPP. These include the choice 
of the atomic configuration, the core radius, the use of 
nonlinear core correction (NLCC) 31 etc. Furthermore, 
since the Slater average might change the atomic orbital 
energies of the NCPP, we do not need to keep the all 
electron orbital energies as the NCPP energies during its 
construction. Thus, we are free to adjust the original 
atomic orbital energies in NCPP construction in order to 
optimize OEPP. We choose Ga as an example to show 
the effect of above settings and adjustments on the im¬ 
provement of OEPP accuracy. The testing results are 
shown in Table I, including the eigenvalue of each orbital 
e s and e p , the pseudo-atom total energy E to t, the inte¬ 
grated densities of s and p orbitals inside the core radii 
(p c s and pp, and the differences of integrated total density 
inside core radii for NLPP and OEPP (8 P ). For these val¬ 
ues, the smaller the differences between OEPP and the 
corresponding NCPP, the more accurate the OEPP is. 
As shown in Table I, we find for elements like Ga that 
has very small <5o, OEPP works well. It reproduces the 
total energy and the orbital energies at the constructed 
atomic configurations. The variation of the constructing 
configuration and the other parameters such as cutoff 
radii does not significantly affect its accuracy. 

Furthermore, we examine the transferability of the 
OEPPs by calculating the orbital energies and the total 
energies of atomic configurations other than the one at 



FIG. 2: (Color figures available online) Color map on the feasibility 
of constructing highly transferable local pseudopotentials for the 
elements in the periodic table. The color represents the value of S p 
as shown by the color bar. The values of 5 P for transition metals 
and rare-earth metals are set at 3. 


which the NCPP and OEPP are constructed. The results 
of selected elements including Mg. Ga and Sb are shown 
in Table II. The configurations that pseudopotentials are 
constructed are shown in the first row for each element. 
The total energies are referred to the the ground-state 
energies obtained by the same pseudopotential. The re¬ 
sults are also compared with the bulk-derived local pseu¬ 
dopotentials (BLPS). 42 As shown in Table II, we find 
that both OEPP and BLPS work fairly well for Mg and 
Ga. In contrast, OEPP orbital energies and total ener¬ 
gies show large improvements than BLPS for Sb. Our 
test results reveal that the transferability of OEPPs are 
better than BLPS for the elements later then Ga. 

We extend our study to other elements in the periodic 
table. In general, we find OEPP works well for large 
number of elements including almost all s-block elements 
and many p-block elements. For many of them, for ex¬ 
ample Mg, the accuracy and transferability of OEPP is 
as good as NCPP. The important features of OEPP and 
its corresponding NCPP for a number of elements are 
listed and compared in Table SI. The OEPP features 
are also compared with available BLPS in Table S2. As 
shown in these tables, the performance of OEPP and 
BLPS are in general similar, showing Slater average as 
a favorable method of constructing LPP directly from 
electronic structure of atoms. Furthermore, OEPP out¬ 
performs BLPS in many cases, especially for elements 
later than In. 

An important point elucidated by the OEPP con¬ 
struction and test calculations is that the existence of 
LPP is intrinsic to an element. Despite the use of dif¬ 
ferent parameters, atomic configurations and types of 
NCPPs, the resulting LPPs do not perform well for 
certain elements if the NC condition is not preserved. 
We can use the deviation from the NC condition, i.e. 
8 p = f n |noEPp(0 — ?t-ncpp(t) |d 3 r °f an element as a 
good measure for the possibility of constructing LPP with 
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FIG. 3: (Color figures available online) The radial wavefunction of 
s and p orbitals, the total density and the potential profile for Mg 
(left panels) and As (right panels), obtained using the Trouliers- 
Martins pseudopotentials (TM), BLPS and OEPP. 


high transferability. This condition is weaker than the 
NC condition, which require the preservation of charge 
density inside the cutoff radius for each angular momen¬ 
tum channel. Therefore this condition is necessary but 
not sufficient. We also find S p is not very sensitive to 
different construction of NCPPs. We present them as a 
color map in the periodic table (Fig. 2). For those ele¬ 
ments colored in purple, the LPPs with high transferabil¬ 
ity might be constructed. The elements for which LPPs 
are difficult to construct include second row elements, all 
transition metals, f-block metals and some late p-block 
elements. S p values for transition metals and rare earth 
metals actually go beyond the largest value in the range 
and are therefore fixed at 3. 

An immediate question is why OEPP can work for 
many elements, considering the large change of the po¬ 
tentials in each angular momentum channel after mak¬ 
ing the Slater average. We present the wavcfunctions, 
the total densities and the potentials obtained by using 
NCPP and LPP in Fig. 3 for Mg and As. It shows 
clearly that despite the large variation of the potential, 
the wavefunctions of NCPP and OEPP are almost iden¬ 
tical. This is the direct result of the fact that the pseudo- 
wavefunctions are constructed under many constrains, in¬ 
cluding nodeless, identical to full electron wavefunction 
out of core region and the norm-conserving condition. 
While comparing with the BLPS, we find a large differ¬ 
ence between the two LPPs especially for the region that 



FIG. 4: (Color figures available online.) Band structure of 
hep Mg obtained by use of TM-NCPP(Red), OEPP(Green) and 
BLPS (Blue) respectively. 



FIG. 5: (Color figures available online.) Band structure of GaAs 
in ZB structure obtained by use of TM-NCPP(Red), OEPP (Green) 
and BLPS(Blue) respectively. 


is close to the nucleus. 


B. Bulk properties 

We now test the accuracy of OEPPs in real materials. 
Using Mg and GaAs as examples, we first calculate the 
band structures. Despite the very different way of con¬ 
struction, the OEPP and BLPS yield very similar band 
structures for both Mg (Fig. 4) and GaAs (Fig. 5). The 
results also compare well with the NCPP band struc¬ 
tures, indicating the good accuracy of both LPPs. The 
OEPP bands are closer to BLPS than to NCPP results, 
although both local pseudopotentials are constructed by 
totally different methods. 

Next, we test the bulk properties of selected elemen¬ 
tal solids and binary compounds. For each material, we 
first optimize its geometry at 0 GPa. After that, the vol- 
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umes of the cells are changed from 0.95 Vo to 1.05 Vo, 
and the calculated total energies at each volume is fitted 
by Birch-Murnaghan 3rd order EOS 43 to yield the bulk 
modulus. The energies as function of volumes are shown 
in Figs. 6 and 7 for Mg in hexagonal-close-packing (HCP) 
structure and GaAs in Zinc Blende (ZB) structure. The £ 

calculated cohesive energies for many elemental solids or tS 

the enthalpy of formations for a number of binary com- ^ 

pounds are shown in Tables III and IV, together with >> 

the equilibrium volume and the bulk moduli and their o 

pressure derivatives. 



FIG. 6: (Color figures available online.) The equation of states 
(EOS) of hep Mg, using KS-DFT with LDA exchange functional 
and three different pseudopotentials, including Trouliers-Martins 
pseudopotential TM-NCPP (square), OEPP (circle) and BLPS 
(triangle). The inset show the EOS of total energies shifted by 
the equilibrium total energy as functions of the atomic volumes 
scaled by the equilibrium atomic volume. 

As shown in Fig. 6 and Fig. 7, both OEPP equilibrium 
volume and its dependence on the pressure deviates from 
the NCPP results. Similar to band structure, the OEPP 
EOS for Mg is close to that of BLPS. However, while the 
volume is scaled by the equilibrium volume (V/Vo), and 
the energy is aligned at the minimum point (E — E$) : the 
scaled EOS calculated from NCPP, OEPP and BLPS are 
very close and almost overlap, for both Mg and GaAs. 

As shown in Table III, the OEPP cohesive energies of 
s block elements, such as Li, Mg, Na etc, are generally 
lower than the NCPP results. In contrast, for p block el¬ 
ements, the OEPP results are higher in comparison with 
the NCPP results. There is no clear trend for the equilib¬ 
rium volume, but the OEPP results are generally lower 
than the NCPP results, and they compares better for 
s block elements than p block elements. The available 
BLPS equilibrium volumes are closer to the NCPP vol¬ 
umes than OEPP; however, the BLPS potentials are con¬ 
structed by fitting the bulk properties. Similar trend can 
be found for bulk modulus. 

For binary compounds, as shown in Table IV, the 



FIG. 7: (Color figures available online.) The equation of states 
(EOS) of ZB GaAs, using KS-DFT with LDA exchange functional 
and three different pseudopotentials, including Trouliers-Martins 
pseudopotential TM-NCPP (square), OEPP (circle) and BLPS 
(triangle). The inset show the EOS of total energies shifted by 
the equilibrium total energy as functions of the atomic volumes 
scaled by the equilibrium atomic volume. 


OEPP results generally compare better with the NCPP 
results than for elemental solids. Furthermore, the avail¬ 
able BLPS results are often worse than OEPP. This 
might be caused by the fact that the BLPS potentials are 
constructed for the elemental solids of which the chem¬ 
ical environment is very different. On the other hand, 
the construction of the OEPP can maximize the trans¬ 
ferability. 

In order to test the OEPP in the OF-DFT scheme, we 
also performed OFDFT calculations for Mg and A1 crys¬ 
tals. The bulk properties of Mg and A1 in four different 
structures, including simple cubic (SC), body-centered- 
cubic (BCC), face-centered-cubic (FCC) and hexagonal- 
closed-packed (HCP) structures are calculated using both 
LDA and PBE-GGA exchange correlation functionals 
and presented in Table V. As shown in the table, the 
difference between KS and OF-DFT is less significant 
than that between the OEPP and NCPP. However, only 
local pseudopotentials such as OEPP can be used in 
the OF-DFT calculations. Although the absolute val¬ 
ues of OEPP show considerable discrepancy comparing 
with NCPP results, both KS-DFT and OF-DFT using 
OEPP reproduce the correct order for equilibrium vol¬ 
ume, cohesive energy and bulk modulus of Mg and A1 
in different structures. The errors caused by use of local 
pseudopotentials are systematic. This indicates that OF- 
DFT using OEPP can be used for large scale simulation 
involving large number of structural configurations. 










TABLE III: KS-DFT+OEPP predictions of the bulk properties for a number of elements, including the equilibrium volume 
(Vo), the bulk modulus (Bo), the equilibrium total energy (Eo) and the cohesive energy E c . LDA exchange-correlation function 
is used. 


System 

Structures 

PPs 

Eo(eV) 

E c (eV ) 

V 0 (A 3 ) 

V exp (A 3 ) a 

B 0 (GPa) 

Bo 

Li 

bcc 

NCPP 

-8.5265 

-2.0599 

18.5793 


15.0 

3.2992 



OEPP 

-9.2786 

-1.5259 

18.0179 


16.5 

3.5459 

Na 

bcc 

NCPP 

-102.2505 

-1.4174 

33.3281 

39.4933 

9.1 

3.3497 



OEPP 

-102.2933 

-1.4479 

32.4215 


9.5 

0.6122 

K 

bcc 

NCPP 

-18.0151 

-1.0333 

64.8243 

75.2843 

4.6 

3.8632 



OEPP 

-18.0093 

-1.1134 

63.0346 


4.9 

3.9196 

Mg 

hep 

NCPP 

-24.6368 

-1.8192 

20.8056 

23.2400 

39.7 

3.8592 



NCPP(PBE) 


-1.5200 

22.6020 


36.2 




OEPP 

-24.5898 

-1.5455 

22.0320 


36.6 

4.1335 



BLPS 



21.175 


38 


Ca 

fee 

NCPP 

-39.0517 

-2.0845 

39.6116 

43.4819 

20.1 

3.4662 



OEPP 

-38.5468 

-1.8520 

45.8677 


19.8 

4.1649 

A1 

fee 

NCPP 

-57.2070 

-4.2334 

15.5407 

16.6013 

84.1 

4.1211 



OEPP 

-56.7969 

-3.6493 

18.3104 


58.4 

5.1277 



BLPS 



15.623 


84 


Si 

CD 

NCPP 

-108.1010 

-6.0072 

19.4620 

20.0210 

96.3 

4.1806 



OEPP 

-107.7690 

-5.0208 

22.9126 


61.7 

4.2380 

P 

orthorhombic 

NCPP 

-180.4319 

-5.9977 

19.0757 

19.0280 

81.0 

4.4198 



OEPP 

-181.3171 

-5.3429 

23.3012 


61.5 

3.8757 



BLPS 



13.957 


133 


Ga 

orthorhombic 

NCPP 

-106.9611 

-3.6216 

18.2695 

19.4690 

66.6 

5.7115 



OEPP 

-107.4471 

-4.0423 

16.1598 


83.4 

4.7387 



BLPS 



17.232 


60 


Ge 

CD 

NCPP 

-108.8499 

-4.9892 

22.9649 

22.6350 

65.0 

4.8198 



OEPP 

-110.0396 

-5.3931 

20.8280 


69.7 

-20.0241 

As 

trigonal 

NCPP 

-255.5401 

-5.2602 

21.3964 

21.5210 

80.8 

4.1800 



OEPP 

-256.3345 

-5.7947 

18.8517 


91.6 

4.2379 



OEPP(PBE) 

-255.7398 

-4.9241 

20.3134 


76.0 

4.1581 



BLPS 



20.033 


77 


Se 

trigonal 

NCPP 

-258.6114 

-4.0206 

25.8296 

27.2610 

63.0 

4.2795 



OEPP 

-259.5605 

-4.3883 

23.6493 


69.8 

4.1601 

Br 

orthorhombic 

NCPP 

-367.0428 

-2.0395 

32.8105 


32.6 

4.5644 



OEPP 

-367.9349 

-2.1755 

31.5890 


34.7 

4.4997 

In 

tetragonal 

NCPP 

-209.6031 

-3.3316 

22.6398 

26.1585 

56.9 

5.2405 



OEPP 

-210.7787 

-4.3875 

16.6604 


94.9 

5.6000 



BLPS 



20.052 


64 


Sb 

trigonal 

NCPP 

-153.6718 

-4.9875 

28.4137 

30.2060 

63.6 

4.4383 



OEPP 

-154.0335 

-5.9311 

23.0326 


81.8 

3.1456 



BLPS 



26.816 


63 


Te 

trigonal 

NCPP 

-224.5969 

-3.8268 

32.0412 

33.9250 

56.0 

4.5338 



OEPP 

-225.0817 

-4.3429 

26.2784 


67.4 

3.0025 

Zn 

hep 

NCPP 

-230.8038 

-2.0237 

12.4261 

15.212 

96.5 

4.7049 



OEPP 

-231.8834 

-2.7594 

10.1611 


129.1 

4.5924 


44 


From Ref. 
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TABLE IV: KS-DFT+OEPP predictions of the bulk properties for selected binary compounds, including the equilibrium volume 
(Vo), the bulk modulus (-Bo), the equilibrium total energy (Bo) and the formation energy Ef. LDA is used for the exchange 
correlation functional. . All the BLPS results ref. to 11 


System 

Structures 

PPs 

Eo(eV) 

E c (eV) 

Vo (A 3 ) 

V exp (A 3 ) a 

B 0 (GPa) 

Bo 

GaP 

CD 

NCPP 

-288.314 

-1.519 

39.228 

40.481 

91.5 

4.468 



OEPP 

-289.427 

-1.442 

39.956 


79.2 

4.472 



BLPS 



37.646 


80 


GaAs 

CD 

NCPP 

-363.246 

-0.744 

44.310 

45.166 

73.5 

4.483 



OEPP 

-364.558 

-0.776 

40.105 


86.7 

4.590 



BLPS 



40.634 


75 

4.472 

GaSb 

CD 

NCPP 

-260.885 

-0.252 

54.523 

56.617 

57.5 

4.592 



OEPP 

-261.496 

-0.157 

48.918 


71.3 

4.576 



BLPS 



52.488 


56 


MgSe 

fee 

NCPP 

-285.941 

-2.693 

38.611 

40.492 

66.8 

3.998 



OEPP 

-286.394 

-2.244 

40.089 


66.6 

4.156 

MgTe 

hep 

NCPP 

-252.944 

-3.710 

62.880 

64.846 

38.8 

4.086 



OEPP 

-251.879 

-2.208 

62.068 


41.9 

4.134 

ZnSe 

CD 

NCPP 

-491.150 

-1.735 

42.077 

45.513 

76.2 




OEPP 

-493.377 

-1.933 

35.270 


99.7 

4.338 

ZnTe 

CD 

NCPP 

-456.335 

-0.934 

52.806 

56.773 

57.6 

4.478 



OEPP 

-458.554 

-1.589 

43.940 


80.2 

4.772 

AlAs 

CD 

NCPP 

-313.728 

-0.981 

44.715 

45.383 

73.8 

4.162 



OEPP 

-313.815 

-0.684 

46.128 


70.5 

4.297 



BLPS 



43.616 


80 



‘From Ref 44 
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V. CONCLUSION 


the norm-conserving condition is 


In summary, we proposed a systematic scheme of con¬ 
structing local pseudopotentials directly from the elec¬ 
tronic structure of atoms. This scheme is based on the 
optimized effective potential method and is found to be 
successful in generating local pseudopotnetials for large 
number of elements, with the accuracy and transferabil¬ 
ity close to the first principles pseudopotentials. For most 
of the elements in s and p-block except the second row, 
the LPP can be constructed and the test on real mate¬ 
rials show that they can yield properties comparable to 
the empirical local pseudopotentials that are subtracted 
from the density functional calculations of the solid. 

The test results for atoms and bulk materials also show 
that once the OEPP works well for the elements, it can 
also work well for the compounds formed by these ele¬ 
ments. For many other elements, OEPP may yield large 
errors. In principle, the bulk properties can be restored 
by adjusting the construction parameters or adding cor¬ 
rections to the OEPP. However, it usually leads to lo¬ 
cal pseudopotentials that only work well at the chemical 
condition that it is fitted. Our practice reveals that the 
existence of a valid local pseudopotential with high trans¬ 
ferability is an intrinsic property of the element. 


Acknowledgments 


[ n FFP (r) d 3 r = [ nf GPP (r)d 3 r. (Al) 
J Cl JCl 

It iw worth to notice the analogy between the NC con¬ 
dition and the OEP method. For OEP, the total density 
is conserved at each point in the real space while in NC 
condition the integration of the density for each orbital 
inside the core sphere is conserved. In both cases, the 
shifts of the orbitals are first order. Denoting the orbitals 
of local and semi-local NC pseudopotentials as ipi (r) and 
ipi (r), and their differences as ipi (r), the NC condition 
can be rewritten as: 


i>t (r) Vi (r) d 3 r - 


(r) ipi (r) d 3 r = 0. (A2) 


Assuming that the changes from semi-local pseudopo¬ 
tential to the local pseudopotential is a perturbation, 
the changes of the wavefunction ipi (r) can be expressed 
through first order perturbation as: 


' 4 >* ( r ) = 


Rewriting NC condition, 


V* O') [v (r') - Vi (r')) Vi (r 1 ) . 3 , 

- tfij (r) d 3 r . 

(A3) 
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as 



(A4) 


ip* (r) £iVi (r) d 3 r + c.c. = 0 


, which can be further transformed to 


(A5) 




0 (A6) 


Appendix A: Norm-conserving condition for local 
pseudopotentials 

Denoting the charge densities of each orbital for local 
and semi-local NC pseudopotentials as n\ FF and nf CPP , 


by employing KS equation 

(-^V 2 + Ui(r)^i(r) =£iVi{r). (A7) 

Using partial integration, we have 


J 


f ip* (r) V 2 ipi (r) d 3 r = ip* (r) Vipi (r) |n - f Vip* (r)Vipi (r) d 3 r (A8) 

Jn Jn 

= ( r ) V Vi (r) |n - Vi/)* (r) <pi (r) \n+ f V*ip* (r) tpi (r) d 3 r. 

Jo. 

The first two terms should be exactly 0 while the NC condition is satisfied because the local pseudopotential will now 
produces wavefunctions that identical to the wavefunctions of all electron potential as well as the NC pseudopotentials 
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TABLE V: Comparing OF-DFT and KS-DFT results of bulk moduli (Bo in GPa), bulk equilibrium volumes (Vo in A 3 ) and 
equilibrium total energies (Eo in eV per atom) for Mg and Al. A number of structure are selected including hep, simple cubit 
(sc), body-centered cubic (bcc) and fee. The energies are refer to the energies of the equilibrium ground state structures, which 
are hep for Mg and ccc for Al. Results obtained by five method combinations are compared, including KS-DFT+NCPP+LDA(or 
PBE ), KS-DFT+OEPP+LDA(or PBE) and OF-DFT+OEPP+LDA (using WGC KEDF). 


Mg 


hep 

SC 

bcc 

fee 

Po 

KS-NCPP (PBE) 

20.835 (22.602) 

24.968 (27.257) 

20.748 (22.560) 

20.970 (22.777) 


KS-OEPP (PBE) 

22.023 (23.606) 

26.391 (28.478) 

22.002 (23.656) 

22.210 (23.818) 


OF-OEPP 

22.225 

26.726 

22.226 

22.333 

a o 

KS-NCPP (PBE) 

3.093 (3.188) 

2.922 (3.010) 

3.462 (3.561) 

4.377 (4.499) 

co 

KS-OEPP (PBE) 

5.027 (5.134) 

3.149 (3.237) 

5.126 (5.206) 

2.977 (3.054) 

3.530 (3.618) 

4.462 (4.569) 


OF-OEPP 

3.154 

5.158 

2.989 

3.542 

4.470 

E 0 

KS-NCPP (PBE) 

0.000 (0.000) 

0.405 (0.389) 

0.026 (0.029) 

0.010 (0.013) 


KS-OEPP (PBE) 

0.000 (0.000) 

0.337 (0.334) 

0.028 (0.030) 

0.007 (0.010) 


OF-OEPP 

0.000 

0.322 

0.024 

0.006 

B 0 

KS-NCPP (PBE) 

41.4 (36.2) 

27.2 (22.6) 

43.6 (35.4) 

38.9 (35.2) 


KS-OEPP (PBE) 

36.5 (33.6) 

22.4 (21.2) 

39.3 (32.7) 

35.9 (32.9) 


OF-OEPP 

35.0 

22.6 

34.4 

34.6 

Al 


fee 

bcc 

SC 

hep 

Vo 

KS-NCPP (PBE) 

15.544 (16.525) 

15.906 (16.945) 

19.055 (20.149) 

15.665 (17.922) 


KS-OEPP (PBE) 

18.029 (19.163) 

18.770 (19.615) 

21.218 (22.324) 

18.488 (20.328) 


OF-OEPP 

18.435 

18.723 

21.683 

18.513 

ao 

KS-NCPP (PBE) 

3.961 (4.044) 

3.168 (3.238) 

2.671 (2.722) 

2.798 (2.858) 

Co 

KS-OEPP (PBE) 

4.162 (4.249) 

3.348 (3.399) 

2.768 (2.817) 

4.618 (4.723) 
2.959 (3.009) 
4.876 (4.927) 


OF-OEPP 

4.192 

3.345 

2.788 

2.960 






4.878 

E 0 

KS-NCPP (PBE) 

0.000 (0.000) 

0.106 (0.094) 

0.401 (0.371) 

0.037 (0.031) 


KS-OEPP (PBE) 

0.000 (0.000) 

0.050 (0.048) 

0.254 (0.246) 

0.014 (0.014) 


OF-OEPP 

0.000 

0.054 

0.223 

0.011 

Bo 

KS-NCPP (PBE) 

83.1 (77.2) 

76.2 (67.3) 

61.3 (55.8) 

79.3 (67.6) 


KS-OEPP (PBE) 

69.4 (54.3) 

54.9 (52.3) 

52.2 (45.8) 

57.6 (50.8) 


OF-OEPP 

67.4 

54.8 

50.6 

56.3 
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outside the core region. Therefore, ip* (r) and Vip* (r) should be 0 at the core sphere and in the region out of the 
core. 

As proved by KLI (see Eqn. 59 in Ref. 25 ), 

- (r) = [v (r) - Vi (r) - (fi - vl)] p* (r) + (e* - v t (r)) ip* (r). (A9) 

Thus, we can transform the NC condition as 


VC ( r ) ( - 2 V +Vi ^ ^ ^ d r 


hi t v ~ 9 V ^ ir) ) ^ ^ + ^ ^ Vi ^ ^ j d r 

= [ {[v(r)~ Vi (r) - (iii - vl)] p* (r) p t (r) 

J n 

+ (ei - Vi (r)) -0* (r) pi (r) + ip* (r) Vi (r) Pi ( r)}d 3 r 

= f {(r) - (r) - (iii - vl)] p* (r) pi (r) d 3 r 

Jn 
= 0 


This means that 


/ rii (r) u (r) dr n, (r) u* (r) dr rii (r) (u* — u*) dr = 0 

J o, J q J n 

Recalling that v\ = f n rii (r) u* (r) dr and u* = f n rii (r) v (r) dr, the above equation can be written as 

(v l - v\) = ( v l - v\) f rn ( r) dr = [v l - v •) fn 
Jn 

Because rb is the partial charge enclosed in the core region and should be smaller than 1, (y l — vf) = 0, 
more clearly shown as: 


(A10) 


(All) 
which can be 


Thus, we prove Eqn. 12 in Section II. 


/ rii (r) (v (r) - Vi (r)) dr = 0. 
Jn 


(A12) 


J 


Appendix B: The electronic structure of (pseudo) atom and related definitions 


For a given set of valence occupancies fi the pseudo valence states are determined by self-consistently solving the 
radial Schrdingier equations associated with the pseudopotential 


1 d 2 1(1 + 1) t/H xc 


2 dr 2 2r 2 


+ 1/ HXC (r) + V^(r) - 


u P s (r) = 0 , 


with the screening potential 


F HXC (r) = F xc [p ps + p{ i ore ; r] + F H [p ps ; r] 


and p ps (r) = fi l' u f S ( r )l • The total energy of the pseudo atom is given by 


where the kinetic energy associated with the pseudo valence states is 


2> ps ] = E /« (V - j™ {^ HXC (0 + K S W| 2 dr) . 


(Bl) 


(B2) 


nOO 

E tot-PS [pPS ] = T [ pP s] + ^xc^ps + -core] + [pPS] + V /z / V) ps (r) |< s (r) | 2 dr , (B3) 

, Jo 


(B4) 
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The above equations refer to Eq.(51-54) in Ref. 30 Here several density related quantities are defined as follows: 


r R c . 


Pi = / K S (ef;0 1 2 dr 


(B5) 


It has the following relation for non-local normal-conserving pseudopotential. 


Pi = / I u i 


PS ( e ? S ; r)\ 2 dr = 


ir(eT e ;r)| 2 dr 


(B6) 


Finally, we define 


S P = / |n PS (r) - »tmW |d 3 r 

Jn 


(B7) 


* Electronic address: mym@calypso.cn 

t Electronic address: miaoms@gmail.com 
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